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Harold C. Euler and Peter E. Wasko 
SUMMARY 


An evaluation of the spherical harmonic fitted geomagnetic field 
equations using Jensen and Whitaker's coefficients and of the inverse 
cube geomagnetic field equation is undertaken in this report. The 
total fields calculated with the former equation are compared with 
those computed with the latter equation and with Vanguard III field 
measurements. 

The total geomagnetic fields were computed at 0, 100, 200, 300, 
400, 600, 800, 1000, 1500, 2000, 3000 and 4000 km and at 1, 2, 4, 8 
and 16 earth radii above a 6371. 2 km spherical earth for colatitudes 
0-° 10', 15°, 30°, 45°, 60°, 75°, and 89° both at 80°W and 100°E longi- 
tudes. For each of these locations the fields were calculated using 

the Jensen and Whitaker's coefficients for S = 1, 2, 3 296, 

where S is the total number of Legendre polynomials arranged in the 
sequence used by Schmidt. The inverse cube fields were also calcu- 
lated for these locations. The Vanguard III field measurements used 
were in the altitude range, 600 - 3800 km, in the geographic area, 
from 2° to 34°N latitude and from 75° to 83° W longitude. 

Various graphs showing the percent truncation levels, S trunca- 
tion levels, and percent deviations of inverse cube fields and measured 
fields from the fields computed for A = 296 Legendre polynomials 
were constructed and interpreted. 

The most important deductions of the evaluation study are sum- 
marized and presented in the section on conclusions. 
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SECTION I. INTRODUCTION 


A determination of the expected geomagnetic field environment 
is important for input to launch vehicle and spacecraft design and 
performance studies. Many control features have torque -producing 
magnetic components which tend to change vehicle spin and orientation. 
The resultant torque on the spacecraft is equal to the cross product of 
the spacecraft magnetic moment vector with the total geomagnetic 
field intensity vector. The geomagnetic field, at the earth's surface, 
is of the order of magnitude of 1/2 Gauss varying with, latitude and 
longitude and decreasing with altitude. The vertical extent of the 
geomagnetosphere, typically about 1 0 5 km, varies with solar activity 
which produces solar flare plasmas and changes in the solar wind 
(Refs. 1 and 2). The geomagnetosphere ranges, in extent, from 
approximately 8 earth radii on the daytime side of the earth to about 
60 earth radii (Ref. 1) on the nighttime side depending upon solar 
activity. 

The inverse cube geomagnetic field equation and a spherical 
harmonic fitted geomagnetic field equation of various orders and 
degrees have been used to compute the fields at various altitudes. 
Jensen and Whitaker produced a spherical harmonic fit to the Epoch 
1955. 0 surface geomagnetic field to degree 24 £.nd order 17. In this 
report an evaluation of these two geomagnetic field equations is under- 
taken by (1) determining the truncation errors due to the omission of 
higher degree and order Legendre polynomials in the Jensen and 
Whitaker Epoch 1955. 0 spherical harmonic fit to the surface geomag- 
netic field and (2) comparing the geomagnetic field values computed 
by the spherical harmonic fitted equation with (a) values computed by 
the inverse cube law, and (b) some values observed by the Vanguard 
III Satellite. 
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SECTION II. GEOMAGNETIC FIELD EQUATIONS 


A. SPHERICAL HARMONIC ANALYSIS OF THE MAIN FIELD 

1. The Geomagnetic Field Equation . The geomagnetic 
potential of internal origin (Ref. 3) is expressed as: 


where: 



m <n 


m = 0 


/ ji\ 
r/ 


n+ 1 


(g n cos mX + 


h n sinm\)P n (cos 9) 

( 1 ) 


X is longitude East 

V is the geomagnetic potential 

0 is the colatitude 

a is the mean radius of the earth 

r is the radial distance from earth's center 

P n (cos 6) are Legendre functions 

n refers to the degree of the Legendre polynomial 

m refers to the order of the Legendre polynomial 


By differentiating Eq. (1), the geomagnetic field components are: 
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where (Ref. 2): X is the northward horizontal component 

Y is the eastward horizontal component 
Z is the downward vertical component 


(2c) 


The total geomagnetic field is accordingly: 
B = \Tx 2 + Y z +~Z 2 


( 3 ) 


2. The Schmidt Normalized Legendre Function in the 
Geomagnetic Potential Equation. The Schmidt normalized Legendre 
functions (Ref. 4) can be derived as follows: 


.m 
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where: m = 0, 1, 2, ...n 

n = 0, 1, 2, ... k Q 

S m = 1, when m = 0 and 6 m = 0, when m ^ 0. 


An example using Eq. (4) to generate the Legendre 
polynomials and Eq. (2) to compute the geomagnetic field component 
contributions, X^ 1 , for n <, 2, and m < n follows: 
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The P n , gj 1 ■ and X n expressions derived from Eqs. (4 - 10) for 
1 < n < 2 and 0 < m < n are shown in Table I. 
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Table I 


An Example of the Schmidt Normalized Polynomials, P™, and the X™ Geomagnetic 

Field Contributions for 1 < n < 2 and 0 < m < n 
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3. The Gaussian Coefficients in the Geomagnetic Potential 

Equation . The Gaussian coefficients, g™ and h^ , in Eq. (1) are 
interpreted as the Schmidt's Gaussian coefficients used by Vestine 
(Ref. 5) following Schmidt's introduction of his normalized Legendre 
functions, P^P. The Schmidt's Gaussian coefficients when multiplied 

li 

by the negative of the normalizing factor C n from Eq. (6) produce 
the Gaussian coefficients as introduced by Gauss (Ref. 5). The 
Jensen and Whitaker's coefficients are similar to the Gaussian 
coefficients used by Gauss. To illustrate the effect of applying the 
normalizing factor to the Jensen and Whitaker's coefficients, these 
coefficients were graphed in Fig. 1 to n = 12 and m = 12 in the 
sequence used by Schmidt. According to this sequence the abscissas 
in Figs. 1-3 acre labeled in Legendre polynomial number. A table is 
inserted in each figure to identify the polynomial number with its 
appropriate degree and order and hence with the corresponding degree 
and order of the Gaussian coefficient and normalizing factor. It can 
be seen in Fig. 1 that the amplitude of the curves tends, to increase 
markedly with increasing n. The corresponding values of C n are 
graphed in Fig. 2 in the same sequence as used in Fig. 1; the corre- 
sponding Gaussian coefficients are shown in Fig. 3. It can be seen in 
Figs. 2 and 3 that the amplitude of the normalizing factor increases 
markedly with n and that the absolute values of the unnormalized 
Gaussian coefficients decrease markedly with increasing n. It will 
be noted later that the contributions to the geomagnetic field by higher 
n- terms also decrease appreciably. 

Of interest is the change in values of the lower degree 
and order unnormalized Gaussian coefficients since Epoch 1835. The 
values of these Gaussian coefficients (Refs. 5 - 10) are plotted in 
Fig. 4. The changes in the g - and h - coefficients are probably due to 
such factors as: (1) the total number of points and geographic area 
used in the curve fit, (2) the number of coefficients used to fit the data, 
(3) a real secular trend, and (4) instrument accuracy. 

4. The Radius -ratio Factor in Equation 2 . The radius -ratio 
a n + 2 

factor — : , where "a" is the mean radius of the earth and "r" is the 

r • 

radial distance from the earth's center through the surface geographic 
point to the altitude desired above the earth's surface, is graphed in 
Fig. 5. It can be seen that this factor decreases rapidly with altitude 
and increasing "n". 
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B. INVERSE CUBE ANALYSIS OF THE MAIN FIELD 


According to the inverse cube geomagnetic field equation, 
the geomagnetic field, B^, at any altitude, h, above the earth's 
surface is 


where "B 0 " is the surface geomagnetic field, "a" is the mean earth's 
radius, and "r" is the radial distance from the earth's center. 


SECTION III. EVALUATION PROCEDURES 


A. The procedure of evaluating the contribution of higher degree 
and order polynomials in Eq. (2) to the totad geomagnetic field in 
Eq. (3) is based on the following percent deviation expression: 



100 % 


where Bj^ is the total geomagnetic field using all of the Jensen and 
Whitaker's coefficients (to n = 24, m = 17) and Bg is the field value at 
some truncation level, S, where n < 24 and m < 17), of the series in 
Eq. (2). The truncation level can be identified either in terms of the 
percent deviation, T, of Bg from Bp, or in terms of the total number, 

S, of Legendre polynomials in the Schmidt's sequence (see table 
inserted in Fig. 1) used to compute Bg. 

B. Similar expressions are used to obtain (1) the percent 
deviation of the inverse cube field value at any one altitude from the 
B^ value at the same altitude starting with zero percent deviation at 
the earth's surface and (2) the percent deviation of the Vanguard III 
geomagnetic field measurement from the Bj^ value at the same location. 
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SECTION IV. PRESENTATION AND DISCUSSION 
OF THE RESULTS 


The total geomagnetic field at 0, 100, 200, 300, 400-, 600, 800, 
1000, 1500, 200.0, 3000, and 4000 km and at 1, 2, 4, 8, and 16 earth 
radii above the earth's surface was computed for the meridian planes, 

80 ° W and 100 °E, for colatitudes 0° 10', 15°, 30°, 45°, 60°, 75°, and 
89° in the northern hemisphere. For these locations the field was 
calculated using the Jensen and Whitaker's coefficients for S = 1 
(simple dipole), 2 (centered dipole), 3, 4, 5 (eccentric dipole), . . . 296, 
where "S" is the total number of Legendre polynomials arranged in the 
Schmidt's sequence. The inverse cube field was also calculated for 
these locations starting with the computed earth's surface value using 
all the Jensen and Whitaker's coefficients (i. e. , for S = 296). 

The total field for the location, 30 °N latitude and 80 °W 
longitude (near to Cape Canaveral, Florida), from the earth's surface 
to 16 earth radii is shown in Figure 6 for S = 1, 2, 3, 4, 5, 10, and 
296 and for the inverse cube relation. It can be seen that: (1) The 

values obtained with the inverse cube relation are larger than those 
using Jensen and Whitaker's coefficients and diverge with altitude, up 
to at least one earth radius, (2) the values at S = 2 and 3 are numerically 
closer to the values at S = 296 than those S = 1, 4, and 5, and (3) the 
values change relatively little from S = 10 to S = 296. 

In Figs. 7-9 the field values are shown as a function of total 
number of Legendre polynomials at various altitudes above the earth's 
surface. It is apparent from the curves in Figs. 7-9 that: (1) The 

values fluctuate to limiting values with increasing S, (2) the magnitude . 
of the fluctuation decreases with increasing S and altitude, and (3) for 
80 °W longitude the values at S = 2 and 3 are numerically closer to the 
limiting values than those for S = 4, 5, 6, and 7 at 60° and 89° 
colatitudes while at 30° colatitude this is not so. The truncation levels 
at 0. 1, 0. 5, 1, 2, and 3 percent are indicated by arrows up to 4000 km. 
The S value for any one percent truncation level, T, is chosen such 
that decreasing S by one will produce a percent value greater than that 
of the truncation level. Generally, the T levels shift to lower S' values 
with increasing altitude. The S values at T = 0. 1, 0. 5, 1, 2, and 3% 
were computed for the altitudes and locations indicated in the first 
sentence of this section and plotted in the appropriate meridional 
cross sections at 80 °W and 100 °E longitudes. S isolines were then 
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drawn. The S-isolines for various percent truncation levels are shown 
in Figs. 10 and 11. From these figures it is seen that: (1) The largest 

S values, S gradients and S variability with latitude at any one T occur 
in the atmospheric surface layer, which varies from about 0 - 2000 km 
at T = 0. 1% to about 0 -400 km at T = 3%, and (2) the S values decrease 
with increasing altitude. 

The percent truncation levels for S = 1, 2, 3, 4, 5, 10, 25, 50, 
100, 200, 242, and 279 were calculated for the altitudes and location 
indicated in the first sentence of this section. The largest T's, 
regardless of latitude, were plotted on altitude versus S cross sections 
for 80 ° W longitude and 100 °E longitude; the T's at the location 30 °N 
latitude, 80 °W longitude, which is near to Cape Canaveral, Florida 
were also plotted on ah altitude versus S cross section. T isolines 
were then drawn and are shown in Fig. 12. The highest T's are 
obtained for low altitudes and S's; the maximum value was T = 28. 1%, 
which was obtained for altitude zero and S = 2 in the 100 °E cross 
section (Fig. 12b). A trough of low T’s at S = 2 and 3, and a ridge at 
S = 4 are present in the 80 8 E cross section. The isolines in Fig. 12 
generally slope to lower S values with increasing altitude. 

The altitude variation of T for various colatitudes at 80 °W 
longitude is shown in Fig. 13 for S = 1, 2, 3, 4, 5, and 10 and of the 
percent deviation for the inverse cube relation in Fig. 14. The T 
values below 4000 km are less for S = 2 and 3, than for S = 1, 4, and 
5 at 45°, 60°, 75°, and 89° colatitudes. The inverse cube percent 
deviation values generally diverge from zero at the earth's surface 
showing positive deviations for colatitudes 45°, 60°, 75°, and 89° and 
on the whole below 4000 km at 80 °W, the geomagnetic field values 
computed with the inverse cube relation are higher than the Bj^ values 
while those computed with a small number of Legendre polynomials are 
lower than the B^ values; the reverse situation occurs for low 
colatitudes, 

A comparison by means of percent deviation of measured total 
geomagnetic fields from the Bj^ values computed for the measurement 
locations was also made, For this purpose, the 1959 Vanguard III 
geomagnetic field data between latitudes 2 S N and 34 B N and longitudes 
75 ®W and 83 ® W were used. The percent deviation data were lumped 
into 200 km layer increments and the 5, 50, and 95 percentile levels 
computed. Five values, which ranged from 1. 63 to 1. 67% In the 



11 


3500 - 3700 km layer were excluded in computing the percentile levels 
since it was felt that these values were not based on representative 
values of the earth's main field. These percentile levels are shown 
in Fig. 15 along with the total range. It can be seen that from 600 to 
3800 km ninenty percent of the computed data lies between -1. 00% and 
0. 30% while all the data lie between -1. 00% and 0. 93%. It appears, 
therefore, that the computed Bl values for latitudes 0 - 45°N and 
longitude 80°W vary from a range of 1. 7 - 2. 7% at 600 km to a range 
of 4. 6 - 9. 1% at 3800 km. In fact, it is not necessary to use the entire 
296 Legendre polynomials to calculate B. The number to be used to 
give as good results should lie between those for T = 2% in Fig. 10b 
and those for T = 1% in Fig. 10c. The maximum number of polynomials 
for T = 1% and T = 2% between 600 and 3800 km from 0 - 45°N are 
respectively 19 and 12. Even the use of the 10 polynomials (see Figs. 

12a and 13f) does not exceed T = 2. 3%. 

By using the data in Fig. 15 as a guide line of the percent deviation 
of measured data from computed data bounding curves of the total geomag 
netic field with altitude were computed. These bounding curves were 
obtained by arbitratily adding 0. 03 Bj^ to the maximum Bj_, value and 
-0.G3BL to the minimum Bl value in order to embrace expected extreme 
variations. The maximum and minimum B]_, values for any one altitude 
were selected from the field values computed at intersection points 
over essentially a 5° latitude -5° longitude worldwide intersection grid. 

In Fig. 16 are shown the curves of the minimum and maximum total 
earth's field with altitude, the Bj_, curve and the inverse cube field 
curve at 30°N, 80°W (near Cape Canaveral, Florida), and the bounding 
or envelope curves. The envelope curves from about 25, 000 km to 
32 earth radii were adjusted to allow for the magnetic field carried by 
a solar flare plasma which, in itself, can reduce or increase the geo- 
magnetic field by approximately 100 gammas (Refs. 11 and 12). 


SECTION V. CONCLUSIONS 

An evaluation of the geomagnetic field equation using Jensen and 
Whitaker's coefficients for Epoch 1955.0 and of the inverse cube 
geomagnetic field equation was undertaken to 16 earth radii above 
the earth's surface for the regions at or near 8CfW and 100°E longi- 
tudes from 0 - 90°N latitudes. From this evaluation it appears that: 
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1. The geomagnetic field equation using the Jensen and 
Whitaker's coefficients to n = 24, m = 17 (568 Gaussian coefficients, 

296 Legendre polynomials) gives results much closer to the Vanguard 
III values (2 - 34 °N, 75-83 0 W) measured from 600 to 3800 km than 
does the inverse cube field equation. The deviation of ±1% (Fig. 15), 
from observed fields, of fields computed with the former equation 
corresponds to a maximum total of 19 Legendre polynomials necessary 
to compute at the truncation level of ±1% (Fig. 10c); the inverse cube 
shows a deviation of approximately 2 to 9% (Fig. 14) from the values 
computed by the former equation in this altitude and geographic region. 

2. Reasonable envelope curves (Fig. 16) to contain the world- 
wide total geomagnetic field variation can be arrived at by (a) adding 
+ 0. 03 Bl> where Bp is the value computed for S = 296, to the world- 
wide maximum Bp j and -0. 03 Bj^ to the worldwide minimum Bp and 
(b) considering the magnetic field carried by the solar flare plasma. 

3. The total field near Cape Canaveral, Florida, (Fig. 6) 
computed with the inverse cube relation is larger than that derived 
from the spherical harmonic fitted equation using Jensen and Whitaker's 
coefficients, and the range between them increases with altitude up to 
at least 1 earth radius. The values for the total number of Legendre 
polynomials, S = 2 and 3, are numerically closer to the values at 

S = 296 than those at S = 1, 4, and 5, and the values change relatively 
little from S = 10 to S = 296. 

4. The total field values as a function of the total number of 
Legendre polynomials S (Figs. 7-9) (a) fluctuate to limiting values 
with increasing S, (b) reduce in amplitude fluctuation with increasing 
S and altitude, and (c) are numerically closer to S = 296 values for 

S = 2 and 3 than for S = 4, 5, 6, and 7 at 1 °N, 30 *N, and 80 "W. 

5. The largest S values, S gradients, and S variability with 
latitude (Figs. 10 and 11) at any one truncation level T occur in the 
atmospheric surface layer, which varies from about 0 - 2000 km at 

T = 0. 1% to about 0 - 400 km at T = 3%, and the S values decrease with 
altitude. 

6. The highest percent truncation levels (T's in Fig. 12) are 
obtained for low altitudes and S's, a trough of low T's at S = 2 and 3 
and a ridge at S = 4 are present in the 80 0 W cross sections and 
absent in the 100 °E cross section, and the T isolines generally slope 
toward lower S values with increasing altitude. 


7. In general, for low latitudes below 4000 km at 80° W 
(Figs. 13 and 14) the geomagnetic field values computed with the 
inverse cube relation are higher than the S = 296 values while those 
computed for small S values are lower than the S = 296 values; the 
reverse situation holds for high latitudes. 
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Fig. 7, a. Total geomagnetic field at various altitudaa above tha earth's surface 
for 30* colatitude, 200* E longitude The arrows and corresponding 
numbers identify the percent truncation levels. 
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Fig. 9. ft. Totftl gaomagnstic fiald at various altitudss abovs tha aarth's aurfaea 
for 60* colatituda. 260* E longituda. Tha arrows and eorrasponding 
numbars idsntify tha pare ant truncation lavals. 
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Cron section of the number of Legendre polynomials necessary to attain the 
truncation level of e0. 1% difference from the total geomagnetic field value 
computed with 296 Legendre polynomials. Jensen and Whitaker's 568 Gaussian 
coefficients (1 £ n S 24; OS m S 17) for Epoch 1955. 0 were used. The isolinee 
are labeled in units of number of Legendre polynomials, 5 
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Fig. 11. a. Cron faction of the number of Legendre polynomials necessary to attain the 
truncation level of *3% difference from the total geomagnetic field value 
computed with 296 Legendre polynomials. Jensen and Whitaker’s 568 Gaussian 
coefficients (l 4 n 5 24; 0 5 m 4 l?) for Epoch 1955.0 were used. The isolines 
are labeled in units of number of Legendre polynomials, S. 
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Fig Il.d. Croat aection of the number of Legendre polynomials necetaary to attain the. 

truncation level of ±0. 5% difference from the total geomagnetic field value 
computed with 296 Legendre polynomial#. Jenaen and Whitaker's 568 Gauaaian 
coefficients (1 5 n £ 24; Oi 17) for Epoch 1955.0 were used. The ieolinea 
are labeled in units of number of Legendre polynomial# , S. 
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Fig. 11. e. Cross section of the number of Legendre polynomials necessary to attain the 
truncation level of *0. 1% difference from the total geomagnetic field value 
computed with 296 Legendre polynomials. Jensen and Whitaker's 568 Gaussian 
coefficients {1 S r S 24; 0 S m i 17} for Epoch 1955.0 were used. The isolines 
are labeled in units of number of Lengendre polynomials, S. 
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Fig. 12.*. Cross section of the maximum percent truncation levels for 
280* E longitude. The isolines are in percent units. 
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Fig. lZ.b. Cross eection of the maximum percent truncation levels for 
100* E longitude. The isolines are in percent units. 
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rig. 12, c. Cross section of the percent truncation levels at 60" colatitude, 

280* E longitude. The isolines are in percent units. 
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Fig. 13. c. 
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Fig. 13. a. Pa re ant Truncation Uvtli at Varioua Colatitude ■ for S ■ 5 
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Fig. 13.f. Percent Truncation Levels at Various Colatitudes 
for S * 10 
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Fig. 14. Percent Deviations at Various Colatitudes of the Geomagnetic Field Computed 
With the Inverse Cube Relation from that Computed for S = 296 
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Fig. 16 a. Mean and Extreme Geomagnetic Field 
Curves for Epoch 1955.0 
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Fig. 16 b. Mein ind Extreme Geomagnetic Field Curvea for Epoch 1955 0 
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